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PLUMEX  II:  A  SECOND  SET  OF  COINCIDENT  RADAR  AND 
ROCKET  OBSERVATIONS  OF  EQUATORIAL  SPREAD-F 


I .  INTRODUCTION 

The  coordinated  measurements  of  equatorial  apread-F 

conducted  during  July  1979  at  the  Kwajalein  Atoll  involved 

the  launch  of  two  instrumented  rocket  payloads  designed  to 

probe  the  detailed  "in  situ"  structure  of  the  turbulent 

ionospheric  plasma.  The  first  launch  operation  (PLUMEX  I; 

17  July  1979;  0031:30.25  LT)  was  conducted  during  the  late 

phase  in  the  development  and  decay  of  spread-F.  The  associated 

results  and  discussions  of  ground-based  and  rocket-borne 
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diagnostics  have  been  presented  in  companion  papers  *  .  The 
second  rocket  (PLUMEX  II:  23  July  1979;  2157:30  LT)  was 
launched  into  the  mid-phase  of  well-developed  spread-F, 
i.e.,  a  ground-based  ionosonde  showed  full  frequency  and 
range  spread  while  the  Altair  radar  maps  of  meter  sise 
irregularity  contours  displayed  backscatter  plumes  that 
penetrated  to  the  topside  F-layer  and  continued  rising  with 
time.  We  present  here  the  initial  PLUMEX  II  observations  of 
radar  plumes  and  the  "in  situ"  measurements  of  the  rocket- 
borne  plasma  probes  and  mass  spectrometer.  The  results  are 
then  compared  with  PLUMEX  I. 


Mmnwript  mbmittcd  July  26, 1980. 
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II 


EXPERIMENTAL  RESULTS 


A.  Ionospheric  Conditions  and  Radar  Maps 

By  2000  hr  LT  on  Che  night  of  Che  PLUMEX  II  launch, 

ionogram8  showed  Chat  Che  virtual  height  (h'F)  of  Che  F- 

layer  had  risen  at  an  average  rate  of  12  meters/sec  Co  an 

altitude  of  350  km.  AC  Chat  point  verCical  drifting  ceased 

* 

and  shortly  thereafter  full-range  spread-F  was  observed. 

The  virtual  height  remained  constant  until  2130  hr  LT,  when 
upward  drifting  again  commenced  at  an  average  rate  of  18 
meters/sec.  With  full-range  spread-F  still  in  effect  and 
with  the  F-layer  still  drifting  at  an  upward  rate  near  18 
meters/sec  the  PLUMEX  II  rocket  was  launched  (0957:30  UT  on 
day  205;  2157:30,  23  July  1979,  LT) . 

Operating  at  155.5  MHz  (radar  backscatter  from  1  m 
field-aligned  Irregularities)  the  Altair  radar  executed 
consecutive  magnetic  east-west  scans  in  a  plane  that  included 
the  upleg  penetration  of  the  rocket  trajectory.  Figure  1 

3 

presents  the  contours  of  constant  backscatter  strength 
plotted  in  increments  of  10  dB,  ranging  from  0  dB  (domain  of 
small  dots)  to  50  dB  (solid  black).  The  first  panel  in 

Figure  1  shows  backscatter  returns  extending  from  350  to 

# 

1000  km,  with  the  most  Intense  region  (50  dB  level)  centered 
at  635  km.  A  measure  of  horizontal  and  vertical  drift 
velocity  can  be  achieved  by  focussing  on  features  "a"  (the 
50  dB  region  at  635  km)  and  "0"  (the  20  dB  region  at  an 


altitude  of  550  km,  and  125  km  veat  of  Altair  in  the  first 
panel) .  In  the  135  sec  time  lapse  between  the  radar  measure¬ 
ments  in  panels  1  and  2,  the  observations  showed  that  (VQ  (hor), 
Va  (vert))  (290  m/sec,  170  m/sec)  and  (V^  (hor),  (vert)) 

^  (300  m/sec,  150  m/sec).  The  results  can  be  interpreted  as 
an  average  magnetic  west-to-east  plume  drift  velocity  of  295 
m/sec  (about  130  a/sec  faster  than  corresponding  observations 
in  PLUMEX  I)  with  a  corresponding  average  upward  vertical 

drift  velocity  equal  to  160  m/sec.  (Note  that  ionograms 
showed  h'F  moving  upward  at  an  average  rate  of  18  m/sec 

throughout  most  of  the  time  encompassed  by  the  radar  maps  in 
Figure  1.  Therefore,  one  could  conclude  that  the  plume  was 
drifting  upward  at  a  rate  of  142  m/sec  relative  to  the  F- 
layer  bottomside.)  In  the  124  second  interval  between 
panels  2  and  3  the  average  magnetic  west-to-east  drift 
velocity  was  considerably  lower,  with  a  rate  at  100  m/sec. 

During  this  same  Interval  the  average  vertical  drift  velocity 
increased  to  270  m/sec  (252  m/sec  relative  to  the  bottomside 
F-region) . 

B.  Rocket  Profile  and  Comparison  with  Radar 
The  rocket  payload  that  was  launched  into  the  spread-F 
conditions  depicted  in  Figure  1  carried,  among  other  instru¬ 
ments,  a  quadrupole  ion  mass  spectrometer  and  a  pair  of 
pulsed  plasma  probes  (See  Ref.  2  for  Instrument  details 
and  the  payload  configuration).  The  pair  of  pulsed 
probes  simultaneously  tracked  ion  and  electron  saturation 

currents  while  generating  conventional  Langmuir  probe 
2 

characteristics  .  Figure  2  displays  the  upleg  measurements 


of  relative  electron  density  es  observed  by  electron  satu¬ 
ration  currents  Ig  and  corroborated  by  simultaneously  measur¬ 
ed  Ion  saturation  currents.  Because  Ion  and  electron  satu¬ 
ration  currents  have  significantly  different  sensitivities 
to  velocity,  sheath  and  magnetic  field  effects,  variations 
in  one  polarity  current  not  corroborated  by  the  other  were 
attributed  to  the  various  aspect  sensitivities  and  excluded 
from  Figure  2.  This  approach  not  only  facilitates  quick-look 
analysis  and  established  credibility  in  the  interpretation 
of  the  curves  as  relative  electron  density  profiles,  but  also 
proved  to  be  a  valuable  technique  for  eliminating  the  other¬ 
wise  degrading  effects  of  uncontrolled  payload  tumple  and  ACS 
jet  firings  that  resulted  from  the  failure  of  payload  sepa¬ 
ration  from  the  rocket  motor  in  PLUMEX  II.  The  ordinate  in 
Figure  2  has  a  linear  scale  for  time-after-launch  with 

altitude  superimposed  at  SO  second  increments. 

The  profile  shows  that  the  payload  entered  the  very 

bottom  of  the  F-layer  at  t  -  145  sec  (Z  340  km)  ,  a  point 
designated  "A"  in  Figure  2.  From  there  to  an  apogee  near 
581  km  (t  -  390  sec)  the  "in  situ"  probe  measurements  revealed 
a  number  of  irregular  structures,  all  contiguous  and  co-' 
located  on  the  bottomslde  gradient  between  points  "B"  and 
"C"  in  Figure  2.  From  the  F-peak  to  the  topside  apogee  of 
581  km  there  were  no  "in  situ"  observations  of  macroscale 
depletions  or  smaller  scale  irregularities  as  reported  in 
PLUMEX  I  (Ref.  1). 

The  payload's  upleg  trajectory  has  been  superimposed  on 


the  radar  maps  in  Figure  1.  With  regions  A  through  E  (and 


their  associated  times  of  observation)  identified  on  the 
panel  best  matched  for  time  coincidence.  A  step-wise  comparison 
of  "in  situ"  observations  (Fig.  2)  with  the  radar  maps  shows 
the  following  correlations: 

Point  "A"  in  Figures  1  and  2  corresponds  to  the  lowest 
position  of  the  bottomslde  F-layer,  a  result  which  co¬ 
locates  the  very  bottom  of  the  F-region  radar  returns  with 
the  very  bottom  of  the  F-region  ledge.  Before  PLUMEX  I  and 
II  (an  Identical  correlation  was  found  in  PLUMEX  I)  there 
were  no  data  to  determine  their  co-location  or  lack  thereof. 

The  transition  from  locations  A  to  B  is  marked  by  a  20 
dB  increase  in  radar  backscatter  energy,  while  region  B  to 
C  is  marked  by  a  relatively  constant  domain  of  backscatter 
energy. 

Observations  at  D  and  E  occur  along  the  eastern  "wall" 
of  the  plume  and  encompass  an  altitude  domain  stretching  from 
the  F-layer  peak  near  500  km  to  a  topside  altitude  near  the 
vehicle  apogee  of  580  km. 

Figure  3  presents  the  up-  and  downing  profiles  of 

relative  electron  density  as  measured  by  "in  situ"  probe 

baseline  electron  saturation  currents.  (The  integrity  of 

the  downleg  profile  was  established  by  the  same  procedure 

utilized  in  Figure  2.)  A  comparison  of  the  profiles  shows 

the  downleg  results  more  representative  of  an  undisturbed  F- 

region,  a  result  consistent  with  the  rocket's  west-to-east 

movement  relative  to  the  radar  plume.  This  relative  motion 

(approximately  160  m/sec)  carried  the  rocket  payload  into 

the  undisturbed  F-reglon  just  to  the  east  of  the  observed 
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r*d*r  plume  (i<e>,  Into  the  white  region  on  the  eestern  edge 
of  the  first  psnel  in  Figure  1). 


A  comment  is  in  order  wlt'h  regard  to  the  up-  and  downleg 
altitudes  that  have  been  superimposed  at  50  second  intervals 
in  Figure  3.  In  particular,  attention  is  directed  to  the 
locations  of  the  very  bottom  of  the  two  F-region  ledge 
observations,  340  km  upleg  and  305  km  downleg.  Because  of 
the  9°  dip  in  the  local  geomagnetic  field  the  downleg  observation 
at  650  sec  corresponds  to  338  km  when  projected  along  the 
magnetic  field  line  to  the  payload  position  at  145  sec  up¬ 
leg,  a  result  showing  virtually  no  net  movement  of  the  ledge 
between  the  up  and  downleg  observations.  (We  note  that  at 
the  time  of  launch  ionograms  showed  h'F  rising  at  an  approximate 
rate  of  18  m/sec.  This  vertical  drift  then  reversed  near 
the  middle  of  the  flight,  establishing  consistency  in  rocket 
and  ionogram  results.) 

With  regard  to  ion  composition,  we  note  that  the  molecular 
ion  NO*  was  the  dominant  positive-charge  component  at  altitudes 
up  to  340  km,  the  very  bottom  of  the  F-region  ledge  ("A"  in 
Figure  2) .  At  all  altitudes  above  340  km  0*  was  observed  to 
be  the  dominant  ion  (See  Ref.  4  for  detailed  ion  profiles). 

Ill.  CONCLUSIONS  AND  COMPARISONS  WITH  PLUMEX  I 

PLUMEX  II  was  the  second  in  a  two-rocket  operation  that 
successfully  executed  space-  and  time-coincident  measure¬ 
ments  of  equatorial  spread-F.  Comparison  of  results  and 
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focus  on  differences  In  prevailing  ionospheric  conditions 
can  be  achieved  by  an  analysis  of  the  upleg  electron  density 
profiles  observed  in  the  two  operations  (PLUMEX  I  and  II). 
These  profiles  are  presented  in  Figure  4  with  specific 
features  and  conditions  of  observation  listed  in  Table  1. 


TABLE  1 

A  COMPARISON  OF  CONDITIONS  AND  OBSERVATIONS 
IN  THE  DNA/PLUMEX  CAMPAIGN 


PLUMEX  I 

PLUMEX  II 

Launch  Time  (LT) 

17  July  1979;  0031:30.25 

23  July  1979;  2157:30 

Spread-/' 

Conditions 

Late  time;  topside  radar 
plume  in  decay  phase 

Mid-phase  of  well-developed 
spread-/;  plumes  penetrating 
to  topside  and  rising  with  time 

F J^nux 

375  Km 

510  Km 

/'-bottom 

240  Km 

340  Km 

nf* 

1.3  (106)  cm"5 

L  (Gradient  scale 
length) 

4.0  Km 

8.6  Km 

T. 

1350  (±250)°K 

TBD 

Dominant  /2  Ion 

0* 

0* 

Number  of  Depletions 
(A/V,/7V,°  >  .60) 

4 

Only  Bottomside  Spread-/* 

Maximum  Depletion 
WN?)mx 

0.90 

=  0.75 

(An  estimate  of  bottomside 
macroscale  structure) 

"In  Situ"  Irregularity 
Strength 

±80%  fluctuations 
on  bottomside  gradient 

Fluctuation  levels  are  much 
less  intense  than  PLUMEX  I 

Ion  Signatures  in  Holes 

AW  ratio 

None 

Plume  Penetration 
by  Rocket 

Western  "wall"  and 
plume  "head1 

Eastern  "wall" 

Spread-F  conditions.  PLUMEX  I  was  conducted  during  the 
late  phase  in  the  development  and  decay  of  spread-F.  Specifically, 
the  rocket  was  launched  at  0031:30.25  (LT) ,  a  full  3  hrs 
after  the  occurrence  of  full  range  spread  and  only  30 
minutes  before  its  disappearance  from  ground-based  ionograms. 

Major  plume  features  were  relatively  stable  with  respect  to 
vertical  drifts,  and  the  most  intense  regions  of  radar 
backscatter  were  beginning  to  decay. 

The  PLUMEX  II  conditions  were  substantially  different, 
with  the  rocket  having  been  launched  into  the  mid-phase  of 
well-developed  spread-F,  i.e.,  a  ground-based  ionosonde 
showed  full-frequency  and  range  spread  while  the  Altair 
radar  maps  of  meter  size  irregularity  contours  displayed 
backscatter  plumes  that  penetrated  to  the  topside  F-layer 
and  continued  rising  with  time.  The  payload  was  launched  1** 
hours  after  the  onset  of  full  range  spread  and  an  estimated 
two  hours  before  decay. 

F-region  profiles  and  "in  situ"  irregularities.  In  the 

PLUMEX  I  operation,  a  number  of  major  depletions  (AN  /N°  <  0.90) 

e  e 

were  distributed  throughout  the  F-region  with  the  F-peak 

(F  h  )  at  375  Km  and  N  max-  1.3  (106)  cm"3  (±  10%).  The 
4  max  a 

very  bottom  of  the  F-layer  (F-bottom)  was  at  240  km  and  the 

macroscopic  gradient  scale  length  l"1-  (N0)”1  dN°/dz  of 

e  e 

the  bottomside  ledge  was  4.0  Km.  (The  macroscopic  L  was 
calculated  from  a  zero-order  fit  to  the  bottomside  ledge 
between  10  and  10  N®ax.)  The  electron  energy  distribution 
was  characterized  by  T@  -  1350  (-250)°K  with  no  obvious 
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signatures  of  energy  redistribution  in  and  around  the  depletions. 

The  most  intense  "in  situ"  irregularities  (i.e.,  smaller 

scale  structure  imbedded  in  the  larger  scale  depletions) 

occurred  near  275  km  on  the  bottomslde  ledge  (-  80%  fluctuations 

about  an  estimate  N°) .  The  power  spectral  density  in  this 

region  of  intense  irregularities  on  the  bottomslde  was 

-2.5  2 

dominated  by  a  k  ’  power  law  over  the  intermediate  wave¬ 
length  domain  k  =*  2 tt / 1  km  to  k  =  2ir/25  m.  This  result 

supports  the  role  of  the  collisional  Rayleigh-Taylor  instability 

2  5  6 

in  generating  intermediate  wavelength  irregularities  ’  ’ 
during  the  occurrence  of  equatorial  spread-F  (See  References 
1  and  2  for  additional  details.) 

The  PLUMEX  II  data  have  not  received  the  same  level 
of  analysis  already  applied  to  the  first  operation,  a  result 
largely  due  to  the  complications  of  uncontrolled  payload  tumble 
and  ACS  jet  firings  that  arose  from  a  subsystem  failure  to 
separate  the  science  payload  from  the  second  stage  rocket 
motor.  However,  there  are  a  number  of  conclusions  that  can 
be  drawn  at  this  time.  In  the  upleg  density  profile  of  the 
PLUMEX  II  flight  (Fig.  2),  F-region  irregularities  were 
observed  only  on  the  bottomslde  gradient,  with  the  F-peak 
(F_ti „„)  at  510  and  n-8X  ^  6(105)  cm-3.  In  PLUMEX  II  the 
bottomslde  gradient  was  substantially  softer  than  in  PLUMEX  I 
with  F-bottom  at  340  km  and  L  ■  8.6  km.  A  review  of  analog 
records  (uncorrected  for  tumble  and  jet  firings)  indicates 
that  smaller  scale  "in  situ"  irregularities  were  much  less 
intense  than  corresponding  observations  in  PLUMEX  I.  This 
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result  could  support  current  interpretations  (e.g.,  Ref.  7) 
of  east-west  plume  asymmetry  which  suggest  that  the  western 
wall  of  a  plume  (the  PLUMES  1  case)  is  more  unstable  than 
its  eastern  counterpart  (the  PLUMEX  II  case) . 

Ion  composition.  PLUMEX  I  and  II  provided  the  very 
first  vertical  profiles  of  ion  composition  during  the  occurrence 
of  equatorial  spread-F.  The  data  are  particularly  relevant 
to  the  subject  of  plasma  bubble  formation,  vertical  transport 
and  ion  signatures  of  source  domains  for  topside  F-region 
depletions.  Previous  "in  situ"  measurements  of  ion  composition 
within  plasma  depletions  came  only  from  satellite-borne  mass 
spectrometers0  with  results  that  suggested  that  topside 
holes  formed  at  lower  altitudes,  where  [N0+] ,  to*]  >  [0+] , 
and  rise  to  greater  heights  maintaining  (to  first-order) 
their  original  ion  composition.  In  a  final  analysis,  the 
composition  can  be  different  because  of  dependence  on  the 
bubble  formation  process,  the  source  domain  and  the  bubble 
age  as  manifested  by  time-dependent  chemical  processes. 

In  both  PLUMEX  operations  0+  was  the  dominant  F-regio^ 
ion  down  to  the  very  bottom  of  the  F-ledge  (F-bottom).  At 
altitudes  below  F-bottom,  molecular  ions  dominated  with  a 
scale  height  approaching  infinity.  (Detailed  ion  chemical 
results  and  analyses  tb  t  go  beyond  discussions  presented 
here  are  available  in  Reference  4.)  The  connection  of 
this  observation  to  topside  bubbles  and  instability  mechanisms 
is  as  follows: 


The  Rayleigh-Taylor  and  E  x  B  gradient  drift  instabilities, 
which  have  been  proposed  for  bubble  formation,  require  a 
steep  bottomside  gradient.  However,  the  PLUMEX  results  show 
that  the  molecular  ions  are  dominant  only  at  altitudes  below 
the  bottomside  F-region  gradient,  where  scale  heights  are 
very  large  and  inhibit  the  instability  process.  If  the  bubble 
does  indeed  form  on  the  F-region  ledge,  and  if  it  transports 
only  local  ions  to  higher  altitudes,  then  molecular  ions  will 
never  dominate  topside  depletions.  This  apparent  dichotomy 
has  a  number  of  explanations: 

(a)  Relatively  small  depletions  (i.e.,  AH  /N°  £  0.9) 

e  e 

originating  on  the  F-ledge  with  initially  small  horizontal 
extent  can  transport  the  local  ion  composition  to  higher 
altitudes  by  the  Rayleigh-Taylor  process;  but  molecular  ions 
will  not  be  the  dominant  positive  species.  Instead,  the 
{N+J/[0+]  ratio  will  indicate  the  original  source  domain  as 

a  2. 

suggested  by  McClure  et  al.  and  Szuszczewicz  et  al.  ,  and 
discussed  in  detail  by  Narcisi  and  Szuszczewicz*. 

(b)  For  molecular  ions  to  be  dominant  in  a  topside  F- 
region  depletion  it  appears  that  one  of  two  mechanisms  must 
apply:  (i)  An  initially  small  depletion  (i.e.,  ANe/N°  £ 

0.9)  of  large  horizontal  scale  size  can  result  in  much 
higher  depletion  levels  (e.g..  References  8-10)  by  frlnglmg 
fields  that  drawup  the  lower  densities  and  molecular  ions 
that  populate  altitude  regimes  lower  than  the  sight  of  the 
initial  perturbation.  This  mechanism  has  been  studied 

by  Zalesafc  and  Ossakov^  and  appears  in  substantial  agreement 
with  observations. 
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(c)  An  alternate  mechanism  for  molecular  ion  dominance 

in  topside  depletions  has  been  proposed  by  Chin  and  Straus**. 

They  suggest  that  plasma  bubbles  in  the  nighttime  equatorial 

ionosphere  originate  as  wind  driven  waves  at  one  of  the 

highly  variable  density  gradients  below  200  km,  rather  than 

at  the  bottomslde  F-region  ledge.  Once  the  bubble  is  formed, 

with  the  low  densities  and  molecular  ion  dominance  characterized 

by  the  lower  altitude,  it  can  propogate  into  the  bottomslde 

F-region  and  provide  the  initial  perturbation  required  for 

13-15 

the  onset  of  the  Rayleigh-Taylor  mode  . 

Current  analyses  of  PLUMEX  ion  composition  results 
support  conclusion  "a"  as  the  operating  principle  on  the 
nights  of  "in  situ"  investigations.  This  does  not  negate 
"b"  and  "c"  as  candidates  for  other  conditions,  but  it  does 
leave  open  the  question  for  further  experimental  tests. 
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Fig.  1  -  Backseat ter  intensity  contour  naps  (contours  are  in  10 
dB  Increments  with  darkest  domain,  e.g.,a  ,  at  the  50  dB  level) 
with  superposition  of  the  rocket  trajectory.  The  times  identi¬ 
fied  with  observations  "A"  through  "F"  are  included  for  a  more 
complete  representation  of  temporal  correlation  with  the  success 
ive  radar  scans.  Domains  "a"  and  "3"  are  references  for  plume 
dynamics  (see  text) . 
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g.  2  -  Relative  electron  density  profile  as  measures  simultaneously  by  ion 
nd  electron  saturation  probe  currents  collected  on  the  upleg  trajectory. 


Fig.  3  -  Comparison  of  up-  and  downleg  plasma  profiles 
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I .  INTRODUCTION 


Accumulated  information  regarding  equatorial  spread-F 
phenomena  has  pointed  toward  a  definite  causal  relationship 
between  the  large  scale  depletions  (also  referred  to  as 
holes,  bite-outs,  or  bubbles)  and  range-time-intensity 
observations  of  large  ionospheric  domains  with  strong  radar 
backscatter  returns  from  the  much  smaller  (meter  size) 
irregularities  (called  plumes) .  Woodman  and  LaHoz  (1976) 
and  Kelly,  et  al.  (1976) ,  have  suggested  that  a  plume  was 
due  to  a  rising  bubble  leaving  behind  a  wake  of  short  wave¬ 
length  irregularities.  Another  proposal  (Szuszczewicz, 

1978) ,  arising  from  considerations  of  chemistry  and  transport, 
suggested  that  the  radar  returns  originate  across  the  density 
gradients  at  the  boundaries  of  large  scale  depletions.  This 
concept  is  supported  by  the  work  of  Ossakow,  et  al.  (1979) 
where  it  is  inferred  that  a  bubble  rising  through  the  F- 
layer  will  bifurcate  on  its  topside  and  produce  shorter  and 
shorter  wavelength  irregularities,  either  by  a  cascade  or 
two-step  mechanism.  Experimental  evidence  to  verify  this 
position  has  come  from  an  Altair  radar  experiment  (Tsunoda, 
1980a)  which  showed  that  backscatter  maxima  tend  to  occur  at 
altitudes  corresponding  either  to  the  electron  density 
minima  or  the  upper  wall  of  the  plasma  depletion.  More 
recently,  Tsunoda  (1980b)  concluded  that  during  the  decay 
phase  of  meter  scale  backscatter  plumes,  the  radar  returns 
were  maximum  on  the  upper  walls  of  the  plasma  depletions. 

Not*:  Manuscript  submitted  February  22,  1980. 
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Efforts  to  examine  the  exact  relationship  between  radar 
plumes  and  ionospheric  depletions  by  performing  simultaneous 
"in  situ"  and  ground-based  radar  observations  (Kelly,  et  al. 

1976;  Morse,  et  al.  1977)  have  been  limited  to  conditions  of 
bottomside  spread-F  and  required  extrapolations  in  space  and 
time  to  establish  correlations.  As  expected  for  bottomside 
spread-F,  the  "in  situ”  probes  only  observed  plasma  fluctuations 
along  the  portion  of  the  trajectory  below  the  F-layer  peak 
density  while  the  ionosphere  above  the  peak  was  quite  smooth. 

The  ion  composition  within  the  depletions  has  also  been 
the  subject  of  a  number  of  investigations.  Typically, 
satellite  mass  spectrometric  observations  (Brinton  et  al. 

1975;  McClure  et  al,  1977;  Szuszczewicz,  1978)  have  shown 
that  the  ion  composition  can  be  vastly  different  inside  and 
outside  the  bite-outs.  Fe+  ions  may  be  enhanced  or  depleted, 
with  molecular  ions  usually  more  abundant  inside  the  bite- 
out.  Brinton  et  al.  (1975)  and  McClure  et  al.  (1977)  have 
found  0+  depleted  by  as  much  as  a  factor  of  103  to  a  concentration 
below  that  of  N0+.  The  molecular  ion  NO*  was  found  to  be 
dominant  in  the  0+  depleted  region,  with  the  bite-outs 
varying  from  a  few  kilometers  to  tens  of  kilometers  in 
width.  An  analysis  of  the  Atmospheric  Explorer-C  data 
(Szuszczewicz,  1978)  suggested  that  a  given  chemical  volume 
on  the  bottomside  F-layer  (lNO+J,  [02+J  > [0+] )  could  move 
upward  through  a  stationary  neutral  atmosphere  and 
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appear  at  higher  altitudes  as  a  bite-out  in  the  local  plasma 
density.  As  the  bottomside  F-region  plasma  cell  moved  up¬ 
ward,  the  relative  magnitudes  of  its  ionic  components  would 
depend  on  transit  time  and  on  altitude  through  the  height 
distribution  of  the  neutral  gases.  This  model  was  consistent 
with  the  satellite  observations  as  well  as  the  computational 
work  of  Scannapieco  and  Ossakow  (1976) . 

In  a  continuing  effort  to  understand  the  detailed 
relationships  involving  large  scale  plasma  depletions, 
meter-size  irregularities  and  associated  ion-chemical 
signatures,  a  rocket  payload  instrumented  with  a  plasma 
diagnostics  complement  (plasma  probes,  electric  field 
sensors,  mass  spectrometer  and  a  two- frequency  beacon 
experiment)  was  launched  into  the  topside  F-region  ionosphere 
above  Roi-Namur  in  the  Kwajalein  Atoll  (4.3°  N  dip  latitude). 
The  investigation  was  part  of  a  major  effort  which  coordinated 
rocket  and  Altair  radar  observations  with  bottomside  soundings 
and  ground-based  photometric  measurements  of  F-region  winds. 

We  present  here  the  initial  coordinated  observations  of  the 
radar  plumes  and  the  "in  situ"  measurements  of  the  rocket- 
borne  plasma  probes  and  mass  spectrometer. 

II.  EXPERIMENTAL  RESULTS 
A.  Ionospheric  Conditions  and  Radar  Maps 

By  2100  hr  LT  on  the  night  of  the  rocket  launch, 
ionograms  showed  that  the  nominal  bottomside  of  the  F-layer 
had  risen  to  an  altitude  of  400  km.  At  that  point,  the  F- 
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layer  began  drifting  downward  with  an  almost  immediate 
occurrence  of  spread-F.  The  downward  drifting  continued  (as 
did  the  spread-F)  at  an  approximate  average  velocity  of 
lOm/sec  with  the  bottomside  F-layer  having  descended  to  an 
altitude  near  270  km  when  the  rocket  was  launched  (12:31:30 
UT  on  day  198;  00:31:30,  17  July  1979,  LT) . 

Operating  at  155.5  Mhz  radar  backscatter  returns  from 
1  m  ionospheric  irregularities)  the  Altair  radar  executed 
consecutive  east-west  scans  in  a  plane  that  included  the 
penetration  of  the  rocket's  upleg  trajectory.  Figure  1 
presents  the  contours  of  constant  backscatter  strength 
plotted  in  10  dB  increments.  (For  details  of  the  Altair 
system  see  Tsunoda  et  al . ,  1979).  The  first  panel  in  Figure 
1  shows  a  backscatter  plume  just  moving  out  of  the  radar's 
eastern-most  field  of  view.  That  plume,  with  its  highest 
and  most  intense  backscatter  region  near  510  km,  is  connected 
to  backscatter  domains  extending  down  to  the  bottomside  of 
the  F-layer.  The  second  panel,  with  a  center  scan  time  137 
seconds  later  than  the  first,  shows  the  intense  backscatter 
region  further  to  the  east,  having  drifted  there  with  an 
approximate  west-east  velocity  of  160m/sec.  In  the  third 
panel  the  plume  has  nearly  moved  completely  out  of  the 
radar's  field  of  view  and  the  intense  region  near  510  km  has 
decayed. 

B.  Rocket  Profile  and  Comparison  with  Radar 

The  rocket  payload  that  was  launched  into  the  spread-F 
conditions  depicted  in  Figure  1  carried  a  quadrupole  ion 
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mass  spectrometer  (from  the  Air  Force  Geophysics  Laboratory) 
a  pair  of  pulsed  plasma  probes  (from  the  Naval  Research 
Laboratory^ ,  vector  electric  field  sensors  (from  Utah  State 
University)  and  a  two-frequency  beacon  experiment  (SRI 
International) .  The  pair  of  pulsed  probes  simultaneously 
tracked  ion  and  electron  saturation  currents  while  generating 
conventional  Langmuir  probe  characteristics  (see  e.g., 

Szuszczewicz  and  Holmes, U 977) .  Figure  2  displays  the  upleg 
measurements  of  relative  electron  density  as  presented  by 
correlated  ion  and  electron  saturation  currents.  The 
ordinate  has  a  linear  scale  for  time-after-launch  with 
altitude  superimposed  at  SO  second  increments.  (Because  ion 
and  electron  saturation  currents  have  significantly  different 
sensitivities  to  velocity,  sheath  and  magnetic  field  effects 
(e.g.,  Szuszczewicz  and  Takacs,  1979)  data  points  not  corroborated 
by  both  polarity  currents  were  attributed  to  the  various 
aspect  sensitivities  and  therefore  were  not  included  in 
Figure  1.  This  approach  facilitates  quick  look  analysis  and 
establishes  credibility  in  the  interpretation  of  the  curves 
as  relative  electron  density  profiles.) 

The  profile  shows  that  the  payload  entered  the  very 
bottom  of  the  F-layer  at  t-  103  sec  (Z-  240  km) .  From  that 
point,  to  an  apogee  near  590  km,  the  "in  situ"  measurements 
revealed  a  number  of  plasma  depletions  depicted  in  the 
figure  as  regions  C,  D-E,  F-G,  H-I,  and  J-K.  The  largest 
depletion  was  in  region  H-I  where  ANe/N°  -  0.85  with  a  half- 


minimum  vertical  extent  approximatley  equal  to  23  km. 

In  the  regions  of  the  large-scale  depletions,  the  "in 
situ”  measurements  also  revealed  much  smaller  scale  irregu¬ 
larities.  The  central  plot  of  "irregularity  intensity"  in 
Figure  2  identifies  the  regions  of  smaller  irregularities 
and  attempts  to  establish  a  preliminary  quantification  for 
their  intensity.  ("Irregularity  intensities”  were  scaled 
directly  from  probe  current  fluctuations  about  an  estimated 
mean.  As  an  illustration,  the  -4.5  to  +4.5  irregularity 
intensity  within  region  C  represents  a  factor  of  9  in  the 
largest  peak-to-peak  fluctuation  measured  in  that  region. 

(If  vehicle  potential,  plasma  temperature  and  ion  composition 
were  constant  during  the  irregularity  measurements,  then 
I  «  Ne>)  More  quantitative  analyses  along  with  power  spectral 
densities  will  be  determined  for  future  publication.)  The 
results  show  that  the  most  intense  irregularities  occurred 
on  the  bottomside  gradient  (region  C)  with  corresponding 
measurements  at  all  other  altitudes  at  a  much  lower  level . 

We  note  that  the  fluctuations  in  the  largest  depletion 
(region  H-I)  are  smaller  than  those  at  "C”.  The  data  also 
indicate  that  the  more  intense  fluctuations  occur  on  positive 
density  gradients  (C,D,E,  and  I). 

The  payload's  upleg  trajectory  has  been  superimposed  on 
the  radar  maps  in  Figure  1  with  domains  A  through  K  (and 
their  associated  times  of  observation)  identified  on 
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the  panel  best  matched  for  time  coincidence  with  the  radar 
results.  A  step-wise  comparison  of  "in  situ"  irregularity 
observations  (Fig.  2)  with  the  radar  maps  reveals  some 
interesting  correlations: 

Point  "A"  corresponds  to  the  lowest  position  of  the 
bottomside  F-layer,  while  ”B"  is  midway  up  the  steep  bottom- 
side  and  very  near  the  point  of  maximum  positive  density 
gradient.  Region  "C"  is  at  the  boundary  of  the  third  highest 
backscatter  level  (30  dB) ,  and  appears  to  represent  the  mid¬ 
phase  development  of  large  scale  Rayleigh-Taylor  turbulence. 
Observations  at  *D",  "E" ,  ”F"  and  NGN  occur  along  the  western 
"wall"  of  the  plume,  and  encompass  an  altitude  domain 
identified  with  the  F-layer  peak.  Point  "G"  represents  the 
payload's  entry  into  the  large  scale  depletion  centered  near 
240  sec  (490  km)  on  the  upleg  trajectory.  The  payload's 
transit  from  "G"  to  "I"  is  marked  by  a  positive  gradient  in 
backscatter  radar  energy,  with  the  maximum  return  occurring 
on  the  topside  (region  "I"  and  above)  of  the  H-I  depletion. 
Above  the  large  scale  depletion,  observations  "K"  and  "J” 
begin  to  track  the  western  "wall”  of  the  plume  in  the  topside 
F-region. 

C.  Ion  Composition 

0+  was  observed  to  be  the  dominant  ion  component  through¬ 


out  the  entire  F-region.  From  points  of  view  focussed  on 
turbulence  and  transport  the  chemical  constituency  of  two 
regions  are  worthy  of  note: 


In  the  H-I  depletion  on  the  topside  F-layer  the  major 
observed  ion  components  were  [0+]  -  0.998  N  ,  (N+l  -  0.002 
N  and  ( IN0+]  +  [0,+])<  10”*N  .  In  the  adjacent  domains 
there  was  a  different  distribution  of  ions,  i.e.,  outside 
the  depletion  we  found  [0+]  -  0.992  N  ,  [N+]  -  0.007  N  and 

c  c 

( [NO+]  *  [0+])<  2  (10~5)  N  ,  a  distribution  typical  of  the 
«  6 

zero-order  ionosphere  at  those  altitudes. 

The  ion  composition  within  the  H-I  depletion  suggests 
that  it  may  have  originated  at  or  near  the  bottomside  F- 

-f  f\,  + 

region  where  [O  ]  -  [0  ]  H_j*  Such  a  region  exists  at  112 
sec  (Z-262  km)  on  the  upleg  trajectory  where  it  was  observed 
that  [NO+]  and  [Oj]  were  1-2%  of  Nfi  and  the  [0+]/[N+]  ratio 
was  nearly  identical  to  that  observed  in  the  H-I  domain. 

This  points  to  N+  as  a  long-lived  tracer  ion  for  bottomside 
source  regions  of  topside  depletions.  The  fact  that  the 
source  region  levels  of  N0+  and  02+  have  not  been  preserved 
in  the  topside  hole  results  from  their  losses  by  dissociative 
recombination  and  a  simultaneous  loss  in  production  by  ion- 
atom  interchange  and  charge  exchange  reactions  since  [N21 
and  [02]  decrease  markedly  with  altitude.  The  longer  it 
takes  a  bottomside  depletion  to  move  upward  into  the  topside  F- 
layer,  the  more  likely  the  elimination  of  molecular  ion 
signatures  when  [0+]  £  N0  >*[N0+]  +  [O^l ) .  In  the  case  of 
the  H-I  depletion,  a  vertical  transport  time  greater  than 
360  seconds  would  account  for  the  molecular  ion  deficiency. 

(To  arrive  at  this  estimate  we  assumed  an  instaneous  dis¬ 
placement  of  the  bottomside  ion  composition  to  the  H-I 
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altitude  and  calculated  that  in  about  6  minutes  the  molecular 
ions  would  decrease  to  a  concentration  less  than  5  cm”'*.) 

This  time  estimate  suggests  an  upper  limit  of  about  600 
m/sec  for  the  depletion's  average  vertical  drift  velocity,  a 
value  which  is  consistent  with  the  wide  range  in  predicted 
bubble  rise  velocities  (Ossakow  and  Chaturvedi,  1978; 

Ossakow  et  al. ,  1979;  Anderson  and  Haerendel,  1979).  (While 
this  conclusion  is  correct  in  its  own  right,  we  note  that 
Altair  data  prior  to  that  shown  in  Fig.  1  reveal  that  the 
backscatter  plume  was  at  the  nominal  altitude  shown  in 
Figure  1  for  more  than  30  minutes.) 

The  second  region  of  special  note  is  "C"  where  it  was 
observed  that  0+  followed  the  intense  plasma  density 
fluctuations  while  the  molecular  ions  NO+  and  02+  (rep¬ 
resenting  @0.5-1.0%  Ne)  did  not.  Such  a  result  has  a 
possible  explanation  in  an  assumption  that  requires  steady 
state  chemical  equilibrium  in  an  0+  dominant  domain. 

(Molecular  ions  in  region  "C”  can  achieve  equilibrium 
concentrations  in  less  than  10  minutes.)  Under  this  condition, 
molecular  ion  concentrations  are  independent  of  0+  and  vary 
only  with  the  scale  height  of  the  neutral  atmospheric 
constituents  N2  and  02>  The  observations  conform  to  this 
model  with  a  standard  zero-order  atmospheric  distribution, 
suggesting  that  neutral  atmospheric  turbulence  is  not  a 
major  source  for  the  observed  plasma  fluctuations  on  the 
bottomside  F-region. 
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D.  North-South  Extent  of  the  Depletions 

Figure  3  presents  the  up-  and  downleg  profiles  of 
relative  electron  density  as  measured  by  "in  situ"  probe 
electron- saturation  currents.  (The  integrity  of  the  downleg 
profile  was  established  by  the  same  procedure  utilized  in 
Figure  2.)  A  comparison  of  the  profiles  shows  very  good 
agreement  in  the  two  observations  of  plasma  depletions. 
Attention  is  directed  to  the  large  scale  depletion  of  the 
topside  (region  H-I) : 

The  up-  and  downleg  measurements  of  the  large  topside 
hole  were  separated  in  time  by  approximately  340  seconds  and 
in  range  by  112  km.  During  this  340  second  interval  Altair 
radar  measurements  of  plume  movement  showed  an  average 
easterly  drift  at  a  160  m/sec  rate  resulting  in  a  total 
eastward  displacement  of  54  km.  Adding  14  km  to  account  for 
depletion  alignment  along  the  magnetic  meridian  (9°E  of  true 
azimuth)  yields  a  calculated  total  E-W  displacement  of  68  km 
between  the  times  of  the  two  rocket  observations  of  the 
hole.  During  this  time  interval  the  payload  had  an  eastward 
range  velocity  approximately  equal  to  215  m/sec,  resulting 
in  an  east-west  separation  in  the  up-  and  downleg  observations 
of  region  H-I  equal  to  73  km.  From  this  we  can  conclude 
near-perfect  up-  and  downleg  targeting  of  the  hole.  The 
agreement  in  the  two  observations  of  the  H-I  domain  there¬ 
fore  suggests  that  the  depletion  is  aligned  with  the  magnetic 
field  for  at  least  112  km.  (Parallel  arguments  dealing  with 
the  depletion  at  (t«Z)  -  (153s,  350  km)  would  suggest  a 
field  alignment  at  least  as  great  as  163  km,  while  the 
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radar  observations  of  Tsunoda  [1980b]  showed  plasma  bubble 
alignment  can  extend  to  1100  km.) 

III.  COMMENTS  AND  CONCLUSIONS 

Space-  and  time  coincident  measurements  of  equatorial 
spread-F  conducted  during  July  1979  at  the  Kwajalein  Atoll 
have  yielded  the  first  definitive  observations  of  small 
scale  irregularities  (@1  meter  )  and  large  scale  plasma 
depletions  as  measured  independently  throughout  the  F-region 
by  ground-based  radar  and  "in  situ"  plasma  instrumentation. 
Preliminary  analysis  of  the  results  leads  to  the  following 
comments  and  conclusions: 

(a)  Within  a  large-scale,  decay-phase,  topside  F-layer 
depletion  where  "in  situ"  irregularities  were  reduced 
(compared  with  the  depletion's  topside  wall),  the  radar's 
backscatter  energy  was  also  reduced  (compared  with  the  top¬ 
side  wall) .  This  result  suggests  the  co-location  of  maximum 
radar  returns  with  the  upper  regions  of  a  depletion  (its 
topside  wall  and  above)  and  not  with  the  depletion  minimum 
or  bottomside  wall. 

(b)  The  "in  situ"  measurements  established  field 
alignment  of  large  scale  depletions  to  distances  at  least  as 
great  as  163  km.  This  result  supports  the  topside  sounder 
data  of  Dyson  and  Benson  (1968) ,  the  airglow  observations  of 
Weber  et  al. ,  (1978) ,  the  recent  radar  measurements  of 
Tsunoda  (1980  a,b)  and  the  assumption  of  depleted  flux  tubes 
in  the  theoretical  considerations  of  Anderson  and  Haerendel 
(1979) . 
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(c)  Ion  composition  measurements  within  a  topside 

depletion  showed  little  evidence  o£  bottomside  molecular 

tracer  ions  (i.e.,  N0+  and  02+) .  This  result  points  to  the 

requirement  for  rapid  bubble  rise  velocities  and/or  low 

4  -3 

plasma  densities  (N  <<  10  cm  )  within  the  hole  if  the 
bottomside  molecular  ion  composition  is  to  be  maintained  as 
bubbles  drift  upward  through  the  F-region  (Szuszczewicz, 

1978) .  However,  the  measurements  did  reveal  N+  as  a  longer 
lived  tracer  ion.  This  helped  identify  the  lower  altitude 
source  region  to  be  on  the  bottomside  F-layer  gradient. 

(d)  Strong  irregularities  on  the  bottomside  F-region 
gradient  showed  0+  following  large  scale  density  fluctuations 
while  the  molecular  ions  NO*  and  02+  were  relatively  constant. 
Preliminary  analysis  of  this  result  suggests  chemical 
equilibrium  and  eliminates  neutral  atmospheric  turbulence  as 

a  major  source  of  the  bottomside  plasma  irregularities. 

(e)  An  estimate  of  bubble  rise  velocity  was  arrived  at 
by  preliminary  ion  chemical  analysis  of  composition  within 
the  hole.  The  analysis  suggests  an  upper  limit  of  600  m/sec 
for  the  average  vertical  velocity  of  an  85%  depleted  domain 
(85%  on  the  topside,  100%  at  the  F-peak)  as  it  drifted 
upward  from  its  bottomside  source  region  near  260  km  to  the 
topside  F-layer  at  490  km.  This  upper  limit  is  consistent 
with  the  wide  range  in  predicted  values  (Ossakow  and  Chaturvedi, 
1978;  Ossakow  et  al . ,  1979;  Anderson  and  Haerendel,  1979). 
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spectral  density  in  this  region  of  intense  irregularities  on  the  bottomside  was  dominated  by 
a  k'2-6  power  law  over  the  intermediate  wavelength  domain  k  ■  2»/l  km  to  k  ■  2v/2.6  m. 
This  result  supports  the  role  of  the  collisional  Rayleigh-Taylor  instability  in  generating 
intermediate  wavelength  irregularities  during  the  occurrence  of  equatorial  spraad-F. 
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EQUATORIAL  SPREAD  F: 

“IN  SITU"  MEASUREMENTS  OF  ELECTRON  DENSITY 
TEMPERATURE  AND  DENSITY  FLUCTUATION  POWER  SPECTRA 

I.  INTRODUCTION 

The  DNA/PLUMEX  rocket  payloads  launched  into  the  equatorial 
ionosphere  during  the  July  1979  campaign  carried  a  plasma 
diagnostics  complement  that  included  a  quadrupole  ion  mass 
spectrometer,  a  pair  of  pulsed  plasma  probes,  vector  electric 
field  sensors  and  a  four-frequency  beacon.  The  pair  of 
pulsed  plasma  probes  not  only  provided  simultaneous  measure¬ 
ments  of  electron  density  Ng,  temperature  Te>  and  density 
fluctuation  power  spectra  Pn(k),  but  also  provided  the 
capability  for  a  running  measurement  of  relative  variations 
in  mean  ion  mass  <M^>. 

The  measurements  of  N  and  T  form  the  basic  information 

e  e 

on  the  laminar  condition  of  the  ionosphere,  allowing  for  the 
determination  of  the  ionospheric  plasma  response  to  varying 
geophysical  conditions  (solar  and  magnetic  activity,  winds, 
gravity  waves,  etc.)  and  the  detection  of  triggering  mecha¬ 
nisms  (e.g.,  steep  density  gradients  in  Ne>  for  ionospheric 
irregularities . 

The  measurements  of  SN^ (—» P^ (k) )  yield  important  test 
Information  for  signal  channel  models  as  well  as  candidate 
instability  mechanisms  (e.g.,  collisional  drift  modes,  E  X  B 
and  Rayleigh-Taylor)  which  might  be  active  in  the  ionospheric 
plasma. 

In  this  paper  we  describe  the  experimental  technique, 

payload  configuration,  launch  scenarios  and  present  additional 

1  2 

experimental  results  which  complement  companion  papers  *  . 
filrnnrn-ir  w«bs*tad  Aim  IS,  1*80 
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II.  TECHNICAL  OVERVIEW 


3 

Introduction .  The  pulsed  plasma  probe  technique  (P  is 

the  designated  acronym)  is  a  unique  diagnostic  tool  capable 

of  high  spatial  and  temporal  resolution  of  plasma  parameters. 

The  instrument  is  a  Langmuir-type  probe  using  a  special 

electronic  procedure  for  generating  the  current-voltage 
3  4 

characteristic  *  .  The  result  is  greatly  improved  reliability 

and  expanded  versatility  in  Langmuir  probe  measurements.  As 

3 

a  diagnostic  tool,  the  P  technique  reduces  commonly  found 
distortions  in  derived  electron  densities  and  energy  distributio 
functions.  A  unique  feature  of  the  technique  is  its  ability 
to  measure  simultaneously  the  electron  temperature,  density, 
and  the  density  fluctuation  power  spectrum.  Successful 

3 

applications  of  the  P  technique  include  not  only  rocket  but 
also  satellite  and  laboratory  beam-plasma  studies^  of 

turbulent  charged-particle  environments. 

Figure  1  shows  two  types  of  probe  operation.  Figure 
1(A)  depicts  a  linear  sawtooth  sweep  voltage  which  represents 
the  conventional  approach  to  Langmuir  probe  operation  wherein 
some  form  of  continuous  voltage  sweep  is  applied  between 
voltage  limits  V_  and  V+.  Fig.  1(B)  shows  the  pulse-modulated 

3 

sweep  which  has  been  utilized  with  P  .  The  voltage  pulses 
which  follow  the  sawtooth  envelope  generate  the  probe's 
current-voltage  characteristic.  During  the  interpulse 
period,  at  constant  voltage  Vg ,  the  collected  probe  current 
Ig  provides  a  direct  measure  of  variations  in  the  probe- 
plasma  system.  The  pulse  duty-cycle  is  short  so  that  the 


probe  rests  st  its  baseline  potential  for  a  period  much 
longer  than  the  pulse  width.  Fig.  2  shows  pulse  and  base¬ 


line  durations  identified  as  T  and  T„ ,  respectively.  So 

on  d 

that  sweep  voltage  transients  will  not  affect  the  value  of 
probe  current,  the  probe  current  is  sampled  at  the  termination 
of  subinterval  T ^  within  the  sweep  pulse  and  at  the  center 
of  the  baseline  interval  Tfi. 

With  T  much  shorter  than  either  T„  or  the  time  constant 
on  B 

4 

of  the  probe  surface  contamination  layer  ,  the  pulse  procedure 
will  maintain  the  surface  condition  and  associated  voltage 
drop  at  a  more  nearly  constant  level  than  when  using  a 
continuous,  slowly-varying  sweep  voltage.  The  resulting 
current-voltage  characteristic  can  then  be  unfolded  from  the 
plasma  density  fluctuations  (SIg)  so  that  the  electron 
temperature  and  density  are  determined  uniquely.  In  addition, 
the  Ig-values  provide  the  raw  data  from  which  density 

fluctuation  power  spectra  are  determined. 

3 

Payload  configuration  and  P  characteristics.  A  pair 
of  pulsed  probes  were  diametrically  extended  from  the 
forward  end  of  the  rocket  payload  (Fig.  3).  The  sensing 
elements  constructed  from  tungsten  wire,  were  isolated  from 
their  extension  booms  by  coaxial  guard  electrodes  driven  at 
the  same  potential  as  the  probes  themselves.  One  of  the 
probes,  defined  as  the  I-probe,  operated  with  Vg^-lv, 
yielding  net  Ion  baseline  current  Ig.  The  other  probe, 
defined  as  the  E-probe,  operated  with  +  2v,  yielding 
net  electron  baseline  current  I*.  Both  probes  generated 
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complete  current-voltage  characteristics  in  Tg  £  400  msec, 
yielding  absolute  values  of  Ne  and  Tg  at  an  approximate  2.5 
Hz  rate.  Maximum  sampling  occurred  at  2048  Hz,  resulting 
in  0.5  meter  spatial  resolution  for  relative  electron  density 
fluctuations  at  a  vehicle  velocity  of  1  km/sec. 

Probe  electrometers  were  set  to  operate  over  a  dynamic 

-10  -4 

range  extending  from  4(10  )  to  2.5  (10  )  amperes,  with 

automatic  switching  over  8  ranges  maintaining  9  bit  accuracy 
for  all  anticipated  ionospheric  conditions.  The  automatic 
ranging  is  best  illustrated  in  Figure  4,  an  actual  in-flight 
analog  record  of  telemetry  channel  outputs  for  the  probe 
currents  and  applied  voltages,  the  roll  magnetometer  for 
magnetic  aspect  determinations,  and  pitch,  yaw  and  roll 
monitors  on  the  ACS  jets.  The  data  sample  presented  in 
Figure  4  was  collected  when  the  payload  was  55  seconds  into 
flight  while  the  probe  electrometers  were  being  driven 
through  a  load  resistor  for  calibration.  The  probes' 
operation  alternated  between  a  fixed-bias  mode  and  a  pulsed- 
sweep  mode,  with  absolute  currents  determined  by  a  simple 
algorithm  which  coupled  the  switching  0-5  v  TM  signal  on  PCM 
channel  26-1  with  the  sweep  current  range  monitor  on  channel 
27-1.  The  record  format  in  Fig.  4  helped  provide  field- 
estimates  of  density  profiles  without  distortions  of  magnetic 
aspect  sensitivities  and  altitude  control  jets. 

III.  RESULTS 

Density  profiles  and  irregularity  structures.  By 
9  P.M.  (LT)  on  the  night  of  the  first  rocket  launch  (PLUMEX  I) 
the  bottomside  of  the  F-region  had  risen  to  an  approximate 

60 


wt.- 


altitude  of  400  km.  The  F-region  then  began  a  downward 
drift  with  a  simultaneous  onset  of  spread-F.  The  downward 
drifting  and  spread-F  conditions  continued,  and  when  the 
bottomside  F-layer  had  descended  to  an  altitude  below  300 
km,  the  rocket  was  launched  (12:31:30  UT  on  day  198;  00:31:30, 
17  July  1979,  LT)  . 

Figure  5  displays  the  upleg  measurements  of  relative 

electron  density  as  presented  by  correlated  ion-  (I*)  and 

•  o 

electron-  saturation  I*"(V+)  currents.  The  ordinate  has  a 

e 

linear  scale  for  time-after-launch  with  altitude  super¬ 
imposed  at  50  second  increments.  Because  ion  and  electron 

saturation  currents  have  significantly  different  sensitivities 

8 

to  velocity,  sheath  and  magnetic  field  effects  ,  variations 
in  I*  and  I^V*)  not  mutually  corroborated  were  attributed 
to  the  various  aspect  sensitivities  and  excluded  from  Figure 
5.  This  approach  facilitated  analysis,  reduced  computer 
time,  and  established  credibility  in  the  interpretation  of 
the  curves  as  relative  electron  density  profiles. 

The  results  in  Figure  5  show  that  a  number  of  major 
depletions  (ANe/N°  <  0.9)  were  distributed  throughout  the  F- 
reglon.  Each  of  the  large  scale  depletions  (identified 
alphabetically)  has  its  own  distribution  of  irregularities, 
illustrated  in  Fig.  6  by  the  expanded  view  of  regions  C,  D, 

H  and  I.  It  is  clear  that  "C"  is  not  a  single  narrow  bite- 
out  but  a  collection  of  rather  large  irregular  structures 
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extending  over  a  total  altitude  domain  of  about  12  It” . 
(Vehicle  velocity  in  region  C  was  2.4  km/sec.)  To  develop  a 
quantitative  view  of  irregularity  fluctuations  observed  in 
the  F-region,  continuous  linear  detrends  were  executed 
throughout  the  entire  upleg  trajectory.  The  variations 
about  those  linear  detrends  were  then  plotted  in  Figure  S  as 
"Irregularity  Intensity",  with  a  maximum  relative  scale  of 
+4.  A  fluctuation  as  great  as  +  4  approximately  represents 
a  +  80%  fluctuation  about  the  linear  detrend.  (Correlation 
of  these  results  with  macroscale  gradients  and  Altair  back- 
scatter  contours  are  discussed  in  a  companion  paper1. 

Absolute  density  and  temperature.  Absolute  values  of 
electron  density  and  temperature  were  determined  by  conven- 

9 

tional  analyses  of  Langmuir  probe  characteristics  with 
appropriate  care  to  eliminate  perturbing  effects  of  surface 
contamination^,  density  fluctuations^'1®  and  magnetic  field 

g 

effects  .  Analysis  of  approximately  25  characteristics  were 

executed  over  the  F-layer  from  340-560  km.  In  each  case  a 

-3  + 

conversion  coefficient  a  =  N@  [cm  ]/Ie(V  )  was  determined 
so  that  the  Ie(V+)  profile  in  Figure  5  could  be  directly 
scaled  to  absolute  electron  densities.  This  procedure 
yielded  a  ■  (5.5  +  0.5)  101®  electrons  cm  ^A  1. 

The  upleg  profile  has  been  reconstructed  in  Figure  7 
with  relative  and  absolute  electron  density  plotted  as  a 
function  of  altitude.  The  result  shows  the  F-peak  at  375 
km,  with  a  maximum  density  of  1.3  (10  )cm  (+  10%). 
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Analysis  of  the  retarding-f ield  region  of  the  saae  set 
of  current-voltage  characteristics  yielded  Tft  ■  (1350  + 
250)°K,  with  no  obvious  signatures  of  electron  energy 
redistribution  in  and  around  the  depletions. 

Intermediate  wavelength  power  spectra.  The  pulsed 
probe  data  provided  an  excellent  opportunity  for  comparison 
with  the  numerical  simulations^  of  the  collislonal  Rayleigh- 
Taylor  (R-T)  instability  at' intermediate  wavelengths. 
Attention  is  focused  on  the  bottomside  F-layer  gradient  and 
region  C,  which  is  believed  representative  of  the  mid-phase 
development  of  the  R-T  process^.  Typically,  computer 
simulations  employ  several  values  for  the  zero-order  gradient 
scale  length 


and  initialize  the  code  with  some  two-dimensional  perturbation 
superimposed.  In  the  work  of  Keskinen,  et  al.1"*-  L  was 
selected  at  5,  10  and  15  km  and  the  perturbation  took  the 
form  12 


(x,y,t»0)  „  _4  _6 

— s -  (10  4)  Sin  (k  y)  Cos  (k  x)  +  2(10  °)  Sin  (2k 

n°  y  x 

e 

with  k  and  k  being  the  horizontal  and  vertical  wavenumbers, 
x  y 

respectively.  Both  k  and  k  were  set  equal  to  2ff/960  «  in 

*  7 

the  simulation.  In  addition,  the  computetion  assumed  that  L 
was  centered  at  300  km. 


y 


y) 
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Under  actual  conditions  encountered  In  PLUMEX  I  (Fig  7), 
the  bottomside  F-layer  gradient  extended  from  240  to  290  km. 

The  question  of  gradient  scale  length  can  be  studied  in 
Figure  8  where  it  is  shown  that  the  bottomside  gradient 
(encompassed  in  the  105-125  sec  time  frame)  is  not  characterized 
by  a  single  value  of  L.  In  region  "C"  (114s  <  t  <  122s)  L 
is  seen  to  vary  between  2  and  10  km,  whereas  adjacent  domains 
(110s  <  t  <  113s  and  122s  <  t  <  126s)  can  be  characterized 
by  L  *  25  km.  We  would  suggest  that  the  adjacent  domains 
are  representative  of  the  zero-order  gradient  scale  length 
and  that  L  *  25  km  would  be  a  more  appropriate  value  in  the 
numerical  simulation. 

In  any  event,  computer  simulations^  with  L  ■  5,  10  and 
IS  km  showed  that  linearly  unstable  modes  saturate  by  non¬ 
linear  generation  of  vertical  modes.  The  results  yield  one¬ 
dimensional  power  laws  (horizontal  and  vertical)  that  vary 
with  a  spectral  index  (  =  n  in  «  k~n)  between  2.0  and  2.5. 

To  explore  this  result  within  the  context  of  region  "C", 
power  spectral  analyses  were  conducted  over  sliding  intervals 

of  2.4  km.  The  results,  presented  in  Figure  9,  show  that 

-2  5 

the  dominant  behavior  is  k  *  over  the  range  k  ■  2ir/lkm  to 
k  ■  2 ir /2.5m.  The  k  behavior  at  t  ■  116.001  sec  is  a 

result  of  the  very  sharp  density  gradient  (see  region  "C" 

Fig.  6)  encompassed  by  the  domain  of  the  spectral  analysis. 

In  general  we  would  conclude  that  our  results  support 
the  numerical  simulations  of  Keskinen,  et  al.*1.  We  do 


point  out  however  thet  e  epectrel  index  verietion  from  2  to 
2.5  is  e  rather  broad  domain.  Further  testing  of  this 
support  can  be  achieved  with  an  L  *  25  kn  simulation  and  a 
downward  drifting  F-layer  model  that  is  more  in  keeping  with 
the  actual  experimental  conditions.  The  F-layer  time- 
history  can  be  important  since  unstable  modes  appear  to 
require  times  in  excess  of  4,000  seconds  to  saturate... a 
time  during  which  the  F-layer  encountered  in  PLUMEX  1 
drifted  downward  in  excess  of  40  km. 

IV.  SUMMARY 

In  each  of  the  two  rocket  operations  (PLUMEX  I  &  II) 

conducted  at  the  Kwajalein  Atoll  during  July  1979,  the  NRL 

pulsed  plasma  probe  performed  flawlessly.  General  results 

concerning  coordinated  rocket  and  radar  measurements  of 

small  and  large  scale  irregularities  have  been  discussed  in 

1  2 

companion  papers  ’  .  Complementary  results  developed  here 
include : 

(a)  In  PLUMEX  I  a  number  of  major  depletions  (AN^/N^c  90Z) 
were  distributed  throughout  the  F-region,  from  the  bottom- 
side  gradient  centered  near  260  km,  through  the  F-peak,  to  a 
topside  altitude  of  500  km.  The  most  Intense  "in  situ" 
Irregularities  occurred  on  the  bottomside  ledge  where  gradient 
scale  lengths  were  found  to  vary  between  2  and  25  km.  The 
power  spectral  density  in  this  region  of  intense  irregularities 
on  the  bottomside  was  dominated  by  a  k  *  power  law  over 
the  intermediate  wavelength  domain  k  -  2ir/lkm  to  k  -  2ir/2.5m. 

The  experimental  conditions  were  reasonably  matched  to  the 


numerical  simulations  of  Keskinen  et  a  1 .  }  ,  and  a  comparison 
of  the  two  resulted  in  general  agreement. 

(b)  In  PLUMEX  I,  the  F-peak  was  at  375  km,  with  Ne*a*  - 
6  —3 

1.3  (10  )  cm  (+  10%).  The  electron  energy  distribution 
was  characterized  by  T@  -  (1350  +  250)°K  with  no  obvious 
signatures  of  energy  redistribution  in  and  around  the 
depletions . 
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Fig.  8  —  PLUMEX  payload  configuration  and  nominal  trajectory. 
(This  figure  1ms  been  adapted  from  a  Saudis  report.) 
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REGION  D 


RELATIVE  FLUCTUATIONS  RELATIVE  FLUCTUATIONS 

Fig.  6  —  Expanded  views  of  density  fluctuations  observed 
in  regions  C,D,H  and  I  of  Fig.  5 
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Fig.  7  —  Relative  and  absolute  profile  of  electron  density  (PLUMEX  I  upleg) 


on  the  bottomride  gradient  of  the  F-region  layer  shown  in  Figs.  5  and  6 
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s  q  —  Power  aoectUl  analyse*  of  density  fluctuations  in  region  **C”  (Fig.  5) 
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77  -  STREET 
NEEDhAV.,  MA  02194 

C1CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 

01 CY  ATTN  J. M.  AEIN 
01CY  ATTN  HANS  WOLFHARD 
01CY  ATTN  JOEL  BEN6ST0N 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEOFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

INTL  TEL  t  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY.  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 

01CY  ATTN  S.R.  GOLDMAN 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL.  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEHIA 
01 CY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  i  SPACE  CO  INC 
P.O.  BOX  504 
SUNNYVALE.  CA  94088 

01CY  ATTN  DEPT  60-12 

01CY  ATTN  O.R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-10 
C1CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 

OiCY  ATTN  D.  CAUFhAN 


KAMA*.  SCIENCES  CORP 
r.C.  SOX  7463 

COLOT. ADO  SPRINGS.  CO  80933 
01CY  ATTN  T.  MEAGHER 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 
01  CY  ATTN  I.  ROTMUELLER 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL.  MA  01854 

01CY  ATTN  K.  BIBL 
01 CY  ATTN  B.  REINISCH 

N.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOWLE 
01 CY  ATTN  P.  WALDRON 
01CY  ATTN  L.  COUGHLIN 
01CY  ATTN  D.  CLARK 
01 CY  ATTN  J.  DAVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BQLSA  AVENUE 
HUNTINGTON  BEACH.  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  HOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01  CY  ATTN  R.W.  HALPRIN 
01 CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01 CY  ATTN  P.  FISCHER 
01 CY  ATTN  W.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  6UTSCHE 
01CY  ATTN  D.  SAPPENFIELD 
01 CY  ATTN  R.  B06USCH 
01 CY  ATTN  RALPH  KILB 
01 CY  ATTN  R.  HENDRICK 
01CY  ATTN  DAVE  SOWLE 
01 CY  ATTN  F.  FAJEK 
C1CY  ATTN  X.  SCHFI6E 
310  ATTN  CONRAD  L.  LONGKIRE 
CIS'  ATTN  WARREN  A.  SCKLt£T£R 
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I-.ITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  KA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 
01CY  ATTN  H.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CL0VERF1ELD  BLVO. 

SANTA  MONICA,  CA  90404 

01CY  ATTN  E.C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONISPHERE  RESEARCH  LAS 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  OYNAMICS  INC. 

P.O.  BOX  3027 
BELLEVUE,  HA  98009 

01CY  ATTN  E.J.  FRENOUH 

PHYSICAL  DYNAMICS  INC. 

P.O.  BOX  1069 
BERKELEY.  CA  94701 

01CY  ATTN  A.  THOMPSON 

RID  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01 CY  ATTN  FORREST  GILMORE 
01 CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  HILLIAM  B.  NRIGHT  JR. 

01 CY  ATTN  HILLIAM  3.  KARZAS 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CV  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 


RAt.D  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  fcOKICA.  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDR0Z1AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATION,  INC. 

P.O.  BOX  2351 
LAJOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01 CY  ATTN  DANIEL  A.  HAMLIN 
01 CY  ATTN  D.  SACHS 
01  CY  ATTN  E.A.  STRAKER 
01 CY  ATTN  CURTUS  A.  SMITH 
01CY  ATTN  JACK  NCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUOBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

01CY  ATTN  DALE  H.  DAVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

01CY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
01 CY  ATTN  SZ 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,.CA  94025 

01 CY  ATTN  OONARO  NEILSON 
01 CY  ATTN  ALAN  BURNS 
01CY  ATTN  6.  SMITH 
01CY  ATTN  L.L.  COBB 
OICY  OAVIO  A.  JOHNSON 
01 CY  ATTN  WALTER  6.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  N.  B«ON 
OICY  ATTN  R.  LIVINGSTON 
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01CY  ATTN  RAY  L.  IEADA8RAKO 
Q1CY  ATTN  C.  CARPENTER 
OICY  ATTN  C.  PRICE 
C1CY  ATTK  J.  PETERSON 
OICY  ATTK  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 
OICY  ATTN  R.  TSUNOOA 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

UNIVERSITY  OF  TOKYO 
ISAS 

KOMABA.  MEGURO-KU 
TOKYO,  JAPAN 

OICY  ATTN  DR.  K.I.  OYAMA 

MAX- PLANCK- INST ITUT 
FUR  PHYSIK  UND  ASTROPHYSIK 
INSTITUT  FUR  EXTRATERRESTRICHE  PHYSIK 
8046  GARCHING  B.  MUNCHEN,  GERMANY 

OICY  ATTN  PROF.  GERHARD  HAERENDEL 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONOO  BEACH,  CA  90278 
OICY  ATTN  R.K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  OEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON.  MA  01802 

OICY  '  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.W.  CARPENTER 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLY 


PLEASE  DISTRIBUTE  ONE  COPY  (EXCEPT  WHERE  NOTED 
OTHERWISE)  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 

ADVANCED  RESEARCH  PROJECTS  AGENCY  (ARPA) 
STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VA  22217 


COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 

OR.  T.  CZUBA 


CAPT  DONALD  M.  LEVINE 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  R.  MEIER  -  CODE  4141 

OR.  TIMOTHY  COfFEY  •  CODE  4000 

OR.  S.  OSSAKOW  -  CODE  4780 

OR.  J.  GOODMAN  -  CODE  4180 

OR.  E.  SZUSZCZEWICZ  -  CODE  4187  (SO  COPIES) 

DIRECTOR  OF  SPACE  AND  ENVIROMtENTAL  LABORATORY 
NOAA 

BOULDER,  CO  80302 

DR.  A.  GLENN  JEAN 
OR.  G.  W.  ADAMS 
OR.  D.  N.  ANDERSON 
OR.  K.  DAVIES 
OR.  R.  F.  DONNELLY 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AIR  FORCE  BASE.  MA  01731 


OR.  T.  ELKINS 
OR.  W.  SNIDER 
MS.  R.  SA6ALYN 
OR.  J.  M.  FORBES 
OR.  T.  J.  KENESHEA 
OR.  J.  APRONS 
DR.  R.  NARCIS1 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 

DR.  M.  B.  McELROY 
OR.  R.  LINDZEN 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK.  PA  16802 

OR.  J.  S.  NISBET 
OR.  P.  R.  ROW  BAUGH 
OR.  D.  E.  BAR  AN 
DR.  L.  A.  CARPENTER 
OR.  M.  LEE 
DR.  R.  DIVANY 
OR.  P.  BENNETT 
DR.  E.  KLEVANS 

UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 
405  HILL6ARD  AVENUE 
LOS  ANGELES.  CA  90024 

DR.  R.  STENZEL 
DR.  F.  V.  CORONITI 
DR.  C.  KENNEL 
OR.  W.  GEKELMAN 

UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CA  94720 

DR.  K.  HUDSON 


OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 


UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 


U.S.  ARMY  ABEROEEN  RESEARCH  AND  DEVELOPMENT 
CENTER  BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD  21001 

DR.  J.  HEIMERL 


DR.  P.  M.  BANKS 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
OR.  K.  BAKER 
OR.  R.  WILLIAMSON 
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CORNELL  UNIVERSITY 
IThACA,  N.Y.  1485C 

DR.  W.  E.  SWARTZ 

DR.  R.  SUDAN 

DR.  D.  FARLEY 

DR.  M.  KELLEY 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

OR.  S.  J.  BAUER/CODE  600 

OR.  R.  HAR TEL/CODE  621 

DR.  R.  GOLDBERG/CODE  912 

DR.  S.  CHANDRA 

OR.  K.  MAEDO 

OR.  R.  BENSON/ CODE  621 

PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  N.J.  08540 

DR.  F.  PERKINS 
OR.  E.  FRIEMAN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  ORIVE 
ARLINGTON,  VA  22202 

OR.  E.  BAUER 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 

OR.  K.  PAPADOPOULOS 
OR.  E.  on 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH.  PA  15213 

OR.  N.  ZABUSKY 
OR.  M.  B10NDI 

DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 

(12  COPIES  IF  OPEN  PUBLICATION 
OTHERWISE  2  COPIES)  12  CY  ATTN  TC 


UNIVERSITY  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

0-10,  MS-664 

LOS  ALAMOS,  NEW  MEXICO  87545 

DR.  M.  PONGRATZ 
OR.  D.  SIMONS 
OR.  G.  BARASCH 
OR.  L.  DUNCAN 
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Washington,  DC  20350 
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